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Abstract 
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1 Introduction 



Ever since the pioneering scattering experiments by Rutherford in the early 1900s, there 
has been the persistent question of how matter would behave during very hard collisions. 
Two examples are the Landau picture where matter would come to a complete stop before 
exploding due to the enormous pressure built-up or the Bjorken picture where target and 
projectile would pass through each other experiencing only partial deceleration in the 
process. Present day high energy nuclear collision experiments provide not only an arena 
for settling this age-old question but also to see if there is any energy dependence that 
would change the collision scenario from one picture to another. 

Another very good reason to study this is related to the original central goal of these 
experiments, which is to recreate deconfmed matter, the so-called quark-gluon plasma, 
believed to have existed only in the early universe. To confirm its existence in the lab- 
oratory at the up-coming Relativistic Heavy Ion Collider (RHIC) at Brookhaven or the 
Large Hadron Collider (LHC) at CERN, one relies heavily on finding evidence amongst 
the many produced particles. These will be affected by the environment in which they 
are produced. In particular the net baryon number left in the central collision region 
or the amount of baryon stopping affects, for example, photon and dilepton production. 
It is therefore of considerable interest to measure this in experiment. Indeed this has 
been done by NA49 at CERN's SPS [|TJ by measuring the net proton or baryon rapidity 
distribution of the final hadrons. 

Recently it was proposed [0 to measure this using the bremsstrahlung associated with 
slowing down of the baryons during the collisions. This method should be simpler than 
measuring hadron rapidities because of the extreme forward focus of the photons emitted 
from relativistically moving targets and projectiles. The measurement is therefore much 
more localized and only photons instead of different types of hadrons need to be detected. 
From a pragmatic point of view, it is better because most hadron detectors are not able 
to cover the full rapidity range. From a physical viewpoint, it can say something about 
the space-time evolution of the collisions 0, [|. This information will definitely not be 
available from rapidity measurements alone. In addition, it can distinguish between the 
Landau and Bjorken picture described above. 

In view of the two different ways of measuring baryon stopping, one would like a way 
to compare the two and be able to communicate between them. In any case, whenever one 
talks about baryon stopping today, one is usually referring to the shape of the rapidity 
distribution. As far as we are aware, there is no attempt to quantify it in any way so that 
the measurement can be put into a more concrete footing. In this paper, we will make 
just such an attempt. Whenever we refer to baryon stopping, we are referring to stopping 
in rapidity space as is traditionally the case and not in velocity space. It will be shown 
that stopping in rapidity and velocity space are entirely different matters and therefore 
it is very important that it is made clear in which it is being described. There can be no 
confusion between the two. 
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2 Quantifying Baryon Stopping or Transparency 

In order to quantify baryon stopping or its inverse, baryon transparency, we aim to find 
a quantity, called S, that satisfies the following requirements. 

(i) It should be equal to unity if there is complete stopping. That is, all baryons end 
up having y — 0. 

(ii) If there is full transparency, and all final baryons move with the original initial 
rapidity yo, this quantity should be zero. 

(iii) One should be able to define its inverse, baryon transparency, T which is related to 
S by the simple relation S = 1 — T so that it has the opposite value in case (i) and 
(ii). 

(iv) S and T should both be equal to half or approximately so when it is clear that there 
is half stopping and half transparency in rapidity space. For example, S = T = 1/2 
when the rapidity distribution is totally flat. 

(v) For different degrees of stopping, S should range between and 1, signifying trans- 
parent to opaque in that order (or the inverse for T). 

(vi) It should be obtainable from both hadron rapidity data and from bremsstrahlung 
measurements. This requirement is essential for bridging the two types of measure- 
ments. 

A quantity that satisfies all these is 

s=1 (1 -^cos 2 g 1/2 ) /■+<*> v{y)p{y) _ 

2t>oCos6'i/2 J-oo 1 — v (y) cos #1/2 

and therefore 

T = (1 -^ 2 cos 2 fli /2 ) /■+«. ^ v(y) p(y) ^ 



2vq cos 0i/2 J -00 1 — v(y) cos 6*1/2 

Here v = v(yo) is the initial velocity in the center of mass frame which is related to the 
initial rapidity y by the general relation between velocity and rapidity 

v(y)=tanhy. (3) 

The p(y) is proportional to the final baryon rapidity distribution dN/dy. For symmetric 
target and projectile, it is defined by 

/+00 
dyp(y)=2. (4) 
-00 
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Accelerator 


v/i/nucleon 
[GeV] 


1-Vq 




yo 


01/2 

(deg) 


SPS I 


17.3 


5.93x10" 


-3 


2.91 


28.83 


SPS II 


19.4 


4.68x10" 


-3 


3.03 


26.97 


RHIC 


200.0 


4.41x10" 


-5 


5.36 


7.93 


LHC 


1500.0 


7.85x10" 


-7 


7.38 


2.87 



Table 1: The value of 9i/ 2 at the various accelerators. 



To see that specifications (i), (ii) and (iii) are satisfied, one can consider two distributions 

P(y) = 2 S(y) (5) 

and 

p(y) = 5(y - y ) + % + y ) (6) 

which correspond to full stopping and full transparency, respectively. It is easily verified 
that one gets S = 1, T = from the first distribution and S = 0, T = 1 from the second. 
Now for specification (iv), the flat rapidity distribution that respects Eq. (^) would be 

P(y) = 1 /yo for \y\ < y (7) 
= otherwise . 

Using the formula 

f -, -7T ; — = — ^ttt + C °o A In cosh y + ln(l — cos 8 tanhw)} (8) 

J l-costftanhy sin 2 6 sin 2 9 1 v y ' 

and the requirement that 5 = T=l/2fora flat distribution, we end up with the equation 

fn 2n a 2 3/> ( X - cos2 0i/2), A + v cos9 1/2 \ f . 
(2 - v ) cos 1/2 - v cos J 9x/2 — In -r 1 - ) = . (9) 

V0 V 1 - W COS 6»i/2 1 

This equation defines the value of the angle #1/2 which is the angle at which a flat rapidity 
distribution will yield S = T = 1/2 for a given t> or y . At this stage, 9i/ 2 is no more 
than a numerical quantity but its meaning will be explained below. In Table J], some 
values of #1/2 have been computed at the various accelerators. The first case at SPS is 
for Pb+Pb and the second is for S+S collisions. The angle becomes smaller as we go to 
higher energies. 

Our specification (iv) requires that rapidity distributions intuitively half way between 
stopping and transparency in rapidity space should be 0.5. To show that this is indeed 
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Figure 1: Simple test cases of rapidity distribution used in the text with 1 > L > X > 0. 

the case, we use a simple test distribution depicted in Fig. |l|. In this figure, we have two 
blocks at the height of 1/(L — X)yo because of Eq. (g) and symmetric about y = 0. The 
parameters L and A allows us the freedom of a range of distributions. For our purpose 
here, we set L = 1 — A and vary the value of A. For any value of A, we have two blocks 
symmetric about ±y /2, and thus they should all have S ~ T ~ 0.5. We tabulated the 
value of S for a set of values of A at the various accelerators in Table [|. We recover the flat 
distribution when A = and we have two delta functions sitting at ±?/o/2 when A = 1/2. 
We see that (vi) is better satisfied as we go to higher and higher energies. Since this way 
of using bremsstrahlung will only be done at RHIC or at LHC, this is good enough and 
we consider this specification met. In any case, they are all fairly close to 0.5. 

One can introduce another special case which is without any doubt half way between 
opaque and transparent. That is 

P(y) = \ (<% - Vo) + 5(y + y ) + 2 5{y )) . (10) 

This distribution is artificial but is ideal for our purpose here. This distribution describes 
half the baryons from the projectile and half from the target sitting at y — and half of 
them from each initial nucleus traveling with the original yo- This is easily worked out to 
give the exact result S = T = 1/2. The purpose of this last distribution is to show that 
we have a sensible definition in hand. 

3 Difference Between Rapidity and Velocity Space 

When one looks for a numerical definition for stopping or transparency, one encounters 
the question of whether this should be in rapidity or velocity space. If one was working 
within Newtonian mechanics, velocity space would have been the automatic choice. This 
is quite logical since one could easily associate stopping with the slowing down of the 
incoming clusters of nucleons. However, it is also traditional to speak of baryon stopping 
while referring implicitly to the shape of dN/dy in rapidity space. So it is in rapidity 
space that we gave this definition in the previous section. It must be stressed that the 
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A 


5 




SPS I 


SPS II 


RHIC 


LHC 





0.50000 


0.50000 


0.50000 


0.50000 


1/5 


0.49552 


0.49581 


0.49921 


0.49987 


1/4 


0.49354 


0.49394 


0.49881 


0.49980 


1/3 


0.49028 


0.49084 


0.49801 


0.49963 


2/5 


0.48820 


0.48886 


0.49836 


0.49946 


1/2 


0.48689 


0.48759 


0.49689 


0.49931 



Table 2: The value of S for special test cases of p(y) depicted in Fig. [T] with L = 1 — A. 
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Figure 2: (a) What a flat p(y) distribution looks like in velocity space and (b) a flat p(v) 
distribution looks opaque in rapidity space at RHIC. 



difference is huge between the two spaces. For example, the flat distribution given in 
Eq. © at RHIC in rapidity space will appear as in Fig. |2| (a) in velocity space because 

p(y) = - => p(v)= 1 . (ii) 

Vo Vo (1 - v ) 

The halfway point of stopping in rapidity space at RHIC would appear to be much closer 
to transparent in velocity space! On the contrary, a flat distribution in velocity space 

p(v) = - =► P (y) = ^ (12) 

vq t> cosh y 

would be more opaque than transparent when one switches to rapidity. This is shown in 
Fig. Q (b). Using this distribution, one finds S ~ 0.97 or T ~ 0.03. For this reason, we 
have insisted that our definition be in rapidity space, conforming to the convention used 
in the heavy ion collision community. 



6 



4 How To Obtain S and T From Bremsstrahlung 
Measurements 



In Sect. 0, the definition given for baryon stopping S and transparency T depended on 
the distribution p(y) which, for symmetric collisions, is related to the net baryon rapidity 
distribution by dN B _g/dy = A p(y) or net proton rapidity distribution by dN p ^ p /dy ~ 
Z p(y). Here Z and A are the atomic and mass number respectively of the incoming target 
or projectile. This assumes that the net proton rapidity distribution is representative of 
or approximately proportional to the net baryon distribution. This does not seem to be 
too bad an assumption if one examines data from SPS Since one could easily calculate 
the value of S and T for some given data of dN/ dy, the important question would be how 
to meet specification (vi) in Sect. 0. The main difficulty is how to obtain the quantity S 
or T from bremsstrahlung measurements. To solve this problem, we now state that S can 
be related to the intensity distribution of the bremsstrahlung emitted in nuclear collisions 
by 



1 - v 2 cos 2 9 1/2 ( 4tt 2 d 2 I 



v 2 sin 2 #i/2 ^aZ 2 dudVt 



9=9 



(13) 



1/2 



The numerical quantity 9i/ 2 has now been given the physical meaning of the opening 
angle from the beam pipe in which direction the soft photons should be measured. To see 
exactly where this angle lies in relation to other directions, one can work out the angle of 
maximum intensity assuming, for example, the case of full stopping S = 1 or complete 
opacity T = using the formulae given in ref. H. In this case, the cosine of this angle 
fl max is related to the initial velocity v by 

cos6 m&x = (2-v 2 y 1/2 . (14) 

It then works out at RHIC to be 6* max = 0.538° and ^ max = 0.072° at LHC. In view of 
the fact that the intensity distribution falls off with the opening angle 9 from the beam 
pipe, these # max are not too far from those 9\/ 2 in Table [l] at the respective accelerators 
so that there will be sufficient intensity at the 9i/ 2 to enable the photon measurements. 

If the reader has not guessed it already, we will now disclose the physical meaning of 
the somewhat mysterious quantity S or T . The inverse of the prefactor to the rapidity 
integral in Eq. (Q) is in fact proportional to the radiation amplitude for full stopping 
M. So one can now see why the rapidity integral itself will always be less than or equal 
to the inverse of this prefactor (but see the next paragraph concerning the soft photon 
requirement). Because of S and T have an origin in the bremsstrahlung intensity dis- 
tribution, they depend on the (charge) hadron rapidity distribution which automatically 
allows them to bridge the two different methods of determining baryon stopping. So far 
we have not mentioned the contribution from charged mesons to photon emissions. They 
could potentially ruin Eq. ([13]). However, pions are by far the most abundant meson 
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L A S T 




L A S T 


1/5 0.991 0.009 
1/3 0.965 0.035 
1/2 0.876 0.123 
3/4 0.661 0.339 
4/5 0.622 0.377 




1/2 1/3 0.700 0.300 
3/4 1/3 0.417 0.583 
4/5 1/3 0.377 0.622 




3/4 1/2 0.229 0.771 
4/5 1/2 0.198 0.801 
1 1/2 0.124 0.876 


1/3 1/5 0.926 0.074 
1/2 1/5 0.800 0.200 
3/4 1/5 0.540 0.460 
4/5 1/5 0.499 0.501 




4/5 3/4 0.045 0.954 
1 3/4 0.018 0.982 




1 4/5 0.011 0.989 



Table 3: The values of S and T for some more general test cases of p(y) depicted in Fig. [L] 
with various L and A at RHIC energies. 



type and they can be positively as well as negatively charged. Therefore contribution to 
bremsstrahlung from mesons cancel out to a large extent || . 

Although it has been expressed in Eq. fll3]) that the intensity distribution should be 
for low energy photons, in practice a few to tens of MeV should be good enough. The 
reason for the soft photon requirement is to remove all nuclear structural dependence as 
well as any potential interference effects. In ref. it was shown that if there were more 
than one component in the acceleration of the nuclear clusters during the collisions, this 
would result in enhancement in and oscillations of the intensity distribution dl/dudQ 
with to. This is a direct result of the interference between the various components in the 
acceleration. For the purpose of our definition, interference would unfortunately just taint 
any value of S obtained from Eq. (p~3|) . Only in the soft uj limit is it free from this type 
of effect. As seen in Fig. 7 of ref. |4j], the intensity at small u is invariant under this. 

5 Examples 

We will now try some example distributions and work out their S and T values. For the 
test distribution in Fig. p] we vary the two parameters L and A to obtain the stopping 
values for the different cases. The first five entries in Table ^ are for central single-block 
distributions centering around y = 0. As the distribution is widened, S decreases towards 
the flat 0.5 value, as expected. The subsequent groups of entries are for symmetric two- 
block distributions shifting progressively away from the center to either side towards 
\y\ = 1. Thus in each group there is the tendency S — > and T — > 1 which is how a 
sensible definition should behave. 

Admittedly these distributions are only test cases designed to show how the numeri- 
cal definition works. However, more realistic distributions can always be approximately 
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reconstructed from thin blocks (strips) of varying heights so the simple distributions used 
do not affect in any way how the definition meets the specifications stipulated in Sec. |2|. 
It may be that in practice our definition would need to be refined but here we have laid 
the groundwork for a simple but sensible numerical definition for baryon stopping. For 
the actual applications of this to real data, and for other ways of using bremsstrahlung 
from more realistic collision models than those used in Hj, we refer the reader to ||. 
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